Presnel laws at curved dielectric interfaces of microresonators 
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We discuss curvature corrections to Fresnel's laws for the reflection and transmission of light at a 
non-planar refractive-index boundary. The reflection coefficients are obtained from the resonances 
of a dielectric disk within a sequential-reflection model. The Goos-Hanchen effect for curved light 
fronts at a planar interface can be adapted to provide a qualitative and quantitative extension of 
the ray model which explains the observed deviations from Fresnel's laws. 
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The fabrication of lasing microresonators [jT],^| and 
opto-mechanical microdevices || has generated a surge 
of interest in the confinement and propagation of light 
in small dielectric structures. Some understanding has 
been achieved from the ray optics of these systems, com- 
plemented by Fresnel's laws of refraction and reflection 
at the interfaces, e.g., in order to identify and describe 
the relevant resonator modes Fresnel's laws give the 
probability of reflection and refraction of plane electro- 
magnetic waves at planar interfaces of media with differ- 
ent refractive index n. Microresonators, however, often 
are so small that the curvature of their boundary cannot 
be neglected. 

In this paper we investigate, in the limit of large 
wavenumbers, the corrections to the Fresnel coefficients 
which appear due to the curvature of the dielectric in- 
terface. The reflection coefficients are obtained via a 
sequential-reflection model B from the resonance widths 
of a microresonator, which are analytically accessible for 
large wavenumbers. The deviations from Fresnel's laws 
are most noticeable around the critical angle for total 
internal reflection, Xc = arcsin(l/n) (where the refrac- 
tive index of the surrounding medium is set to unity) 
and amount in a systematic reduction of the reflection 
probability. 

The reduction of the reflection probability is conven- 
tionally related to tunneling escape at the curved inter- 
face. In view of the previous success of the ray model, 
which often is desired to be retained for its simplicity, we 
provide an alternative qualitative and quantitative expla- 
nation of the deviations by incorporating into this model 
the Goos-Hanchen effect j5]-|l(|. This effect results from 
the interference of rays in a beam with slight variations 
of the angle of incidence and consists of a shift of the ef- 
fective plane of reflection. At a curved interface, it turns 
out that the reflection probability is then reduced be- 
cause the angle of incidence at the effective interface is 
smaller than at the physical interface. 

There is evidence obtained in the context of quantum- 
mechanical scattering problems [ pd]| that incorporating 
the Goos-Hanchen effect is equivalent to a semiclassical 
approximation. Also the effect has been used in Ref. |T2j ] 
to explain the decreased spacing of resonances observed 



in experiments with dielectric spheres, in terms of an 
effective optical size of the cavity which is larger than its 
physical size. Our work can be seen as complementary to 
this previous study, because we are concerned with the 
resonance widths, not only the resonance energies. 

Curvature corrections to Fresnel coefficients have been 
addressed in the past in a number of works, by applying 
various techniques; see for example Refs. E3-M. The 
works closest in spirit to the present paper are those 
which employ the complex ray method, e.g., to describe 
light rays approaching a disk from outside fl5|,|l6| . 

Although we restrict our discussion to circular inter- 
faces, the results for the reflection coefficients should be 
applicable to microresonators of general shape as soon as 
locally the curvature can be approximated as a constant. 

According to Fresnel's laws, a plane electromagnetic 
wave incident on a planar dielectric interface with an- 
gle of incidence x is reflected with the polarization- 
dependent coefficients 



R™ = 



sin 2 (x - rj) 
sin 2 (x + rj) 



R 



TE 



tan 2 (x - rj) 
tan 2 (x + rj) 



(1) 



where TM (TE) signifies transverse polarization of the 
magnetic (electric) field at the interface and rj = 
arcsin(n sinx) is the direction of the refracted beam (ac- 
cording to Snell's law). 

Let us compare the Fresnel coefficients with the reflec- 
tion coefficients at a curved interface with radius of cur- 
vature r c . Their angular dependence can be conveniently 
obtained from the energies and widths of resonance states 
in a two-dimensional circular disk of radius r c . In this ge- 
ometry the two possible polarization directions decouple 
and angular momentum (quantum number m) is con- 
served. We introduce polar coordinates r and (f> and de- 
note the (complex) wavenumber by k. We will concen- 
trate on the case close to geometric optics Re kr c ^> 1 . 

The resonance states are obtained by matching the 
wave field oc J m (nkr)e im ^ inside the disk (with the Bessel 
function J) at r = r c to the wave field oc Hm\kr)e lm ^ 
outside the disk (with the Hankel function H^), where 
the matching conditions follow from Maxwell's equations: 



J m (nkr c ) H$\kr c ) =nJ' m {nkr c ) H${kr c 



(2) 



1 



for TM polarization, and 

nJ m (nkr c )Hg>\kr c ) = J' m {nkr c )H${kr c ) , (3) 

in the TE case (primes denote derivatives) . Given a com- 
plex solution fc, the angle of incidence is obtained from 
the real part by comparing the angular momentum in the 
ray picture and in the wave picture, 



smx 



11 Re kr c ' 
while the reflection probability 



R = exp(4n Im kr c cos x) 



(4) 



(5) 



follows from the imaginary part of k because it deter- 
mines the escape rate, which in turn can be related to R 
by the sequential-reflection model 

Because the discrete set of resonance energies obtained 
from Eqs. (g) and (||) only is meaningful for the disk, let 
us first derive analytical expressions for the resonance 
width as a function of a continuous resonance energy 
which smoothly interpolate between these solutions. It 
is interesting to note Q that one cannot simply expand 
Eqs. (H), (H) in Imfc when k is not close to an exact solu- 
tion. Moreover, for TE polarization Im k would diverge 
at the Brewster angle if it is calculated by inserting the 
Fresnel coefficient ([!]) into Eq. (|J). In order to achieve a 
more accurate expansion we separate out the problematic 
term and cast Eqs. (0) and (H) into the form 



J' m {nkr c ) 
JffiiyihTc) 



T{kr c 



(6) 



with 



T™{x) 



F™(x) 



-2 -r-TE 



F^(x), (7) 



depending on the polarization. In both cases, J-(kr c ) is a 
slowly varying complex function of its argument, and the 
argument can be taken real because Refc ^> |Imfc|. The 
logarithmic derivative of Bessel functions, however, is a 
rapidly fluctuating function, and its dependence on Im k 
has to be worked out carefully. This can be achieved by 
approximation by tangents E^] , 



tan a cos x, 



J' m (nkr c ) 
J m (nkr c ) 

a = m cot x + m X ~ m 7^ — ^ + i n lm kr c cos x, 



(8) 
(9) 



where x is given as a function of Refc by Eq. (jfy. We 
expanded a linearly in Im fc and neglected terms of or- 
der (Refcr c ) _1 . Equation now can be solved for a, 
without any further approximation. ^From the imagi- 
nary part one deduces 




Re kr_ = 50 



TE 



0.3 







L 
[ 


1 1 1 1 

! 






: & : 

' ^D'" 



1 



0.5 1 0.5 

sin x sin % 

FIG. 1. Resonance widths —lmkr c for a dielectric disk 

with Re kr c — 50 and a refractive index n = 3.29 (upper 

panel) and n = 1.3 (lower panel), for TM and TE polarized 

light, as a function of the angular momentum [parameterized 

according to Eq. (^J)]. The analytical result from Eqs. ([to|), 

( p"l| ) (solid curves) is compared with the exact results from 

Eqs. @, §) (squares). 



Im kr r 



ncosx 



Im arctan ■ 



T 

cosx' 



(10) 



Although this can already be taken as the final result, we 
further may insert the uniform approximation |[17[ 



T 



TE 



in cos r/ 



1 



1 



sin 2 r\ 



K 2 /s(z) 



1 



(11) 



and similarly for J 7 ™ = J^ IE /n 2 , with the mod- 
ified Bessel function K, the angle of refraction 77 
(which is a complex number for x > Xc), and z = 
— iRe kr c cos 3 rj/(3 sin 2 rj). In Fig. [I] we illustrate that 
Eqs. (flO), ( pd| ) agree very well with the exact solutions 
of Eqs~@ and (|), even close to the Brewster angle for 
TE polarization, and interpolates smoothly in between. 

The angular dependence of the reflection coefficients 
can now be obtained by combining Eqs. (||) and (pX|), 
giving the final result 



R 



cos x + iJ~ 



cos x — i3~ 



(12) 



Figures g and [| show the result for the two values 
Refcr c = 50 and Refcr c = 150, respectively. Deviations 
from Fresnel's laws (Q) are most visible around the critical 
angle Xc = arcsin(l/n) where the reflection coefficients 
increase rapidly as the regime of total internal reflec- 
tion is approached. The correction consists not only in 
a broadening of the transition interval, but most notably 
also in a shift of this transition region towards higher an- 
gles of incidence, resulting in a systematic reduction of 
the reflection coefficient. 

The deviations from Fresnel's laws in Figs. || and || 
increase as Re kr c is reduced, that is, the more notice- 
able the curvature of the interface is. On the other 
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Re kr„ = 50, n = 1 .54 




Re kr. = 150, n = 1.54 
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FIG. 2. Reflection coefficients R in a dielectric disk with 
Refcr c = 50. The solid curve is the analytical result (p^[), 
which smoothly interpolates between the exact solutions of 
Eqs. (^|), (|^) with real part close to 50 [squares, translated 
into angular-dependent reflection coefficients by Eqs. (^J) and 
(g)]. The dashed curve is the result of incorporating the 
Goos-Hanchen effect into a ray model (assuming for TE po- 
larization that the shift is the same as for TM polarization, for 
reasons explained in the text). The dotted curve is Fresnel's 
law (§. 



hand, in the zero-wavelength limit Re kr c — > oo of ge- 
ometric optics any interface appears planar, and Fres- 
nel's laws should apply without modification. Indeed 
it can be seen that they are recovered from Eq. Jl^ ) 
when the approximation by tangents is also applied to 
the Hankel functions in T, resulting in J^™ = in^ 1 cos 77, 
JF TE = in cos 77. The deviations close to the critical angle 
are directly related to the break-down of this approxi- 
mation when the argument of the Hankel functions be- 
comes smaller than the index. As we will discuss now, 
the curvature corrections to Fresnel's coefficients can be 
obtained within a minimal extension of the ray picture 
when the Goos-Hanchen effect is taken into account (the 
result obtained is given by the dashed curves in Figs. || 
and ||). 

The Goos-Hanchen effect [| 10 refers to the displace- 
ment of the reflected beam when the incident beam con- 
sists of rays with slight variations of the angle of inci- 
dence, and arises because each ray experiences a slightly 
different phase shift when it is reflected. As is illustrated 
in Fig. ||a, the lateral shift along the interface also can be 
interpreted as resulting from a displacement of the effec- 
tive plane of reflection (a signature of this displacement 




FIG. 3. Same as Fig. |2j but for a wavenumber Re kr c = 150. 
The result based on the Goos-Hanchen effect is almost indis- 
tinguishable from the exact result. 



is the increased Wigner delay time which recently has 
been measured at metallic gratings ]l8[.) 

Since at a planar interface the reflection law is not af- 
fected by parallel displacement and fulfilled for the mean 
angles of incidence, the reflection coefficients are not af- 
fected by the Goos-Hanchen effect — the only consequence 
of the slight variation of angles is that their angular de- 
pendence is smeared out. However, the situation changes 
at a curved interface [jllj as shown in Fig. [I|b. The inter- 
section of the incident and the laterally shifted reflected 
ray defines an effective boundary of radius r' c > r c . We 
now can assume that the ray is specularly reflected at the 
effective boundary, resulting in a smaller effective angle of 
incidence V < %, and evaluate the Fresnel reflection co- 
efficients (|l|) at this smaller angle Since the reflection 
probability is then reduced this qualitatively explains the 
observed deviations from Fresnel's laws in Figs. || and ||. 

For a quantitative comparison we need the distribution 
of angles of incidence P(x)i which can be related to the 
radial width oc r c (nRe kr c )~ 2 / 3 of the caustic in the ra- 
dial wave field J m (nkr). (The same distribution of angles 
will also be used to smooth out the reflection coefficients 
as is appropriate even for the planar interface.) The vari- 
ation of angles of incidence arises from the intrinsic cur- 
vature of the beam wave front and also directly from the 
curvature of the interface — note that both mechanisms 
should contribute equally to the Goos-Hanchen effect. 
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FIG. 4. (a) Goos-Hanchen shift at a planar interface. An 
incident beam containing contributions from plane waves of 
slightly different angle of incidence \ appears to be reflected 
at a position that is shifted a distance zgh away from the 
point of incidence. Alternatively, one can think of the beam 
to be reflected at a shifted interface indicated by the dashed 
lines, (b) Goos-Hanchen effect at a curved interface. The 
reflection seems to occur at an interface of larger curvature 
radius r' c > r c under a smaller angle \' < x of incidence. 



We found that for our purposes P(x) is sufficiently well 
approximated by that of a Gaussian beam of half width 
a = r c (nRe fcr c )~ 2 / 3 <C r c . It is good to observe that this 
Gaussian beam geometry does not put us into conflict 
with the finite disk size while at the same time leaving us 
in the paraxial regime nReka = (nRekrc) 1 ^ 3> 1. For 
TM polarization the effective radius of curvature r' c can 
then be calculated by applying the formulas of Ref. |t]] 
(which arc lengthy expressions and hence not given here) . 
The result of this Goos-Hanchcn-cffcct-based approach 
is presented as the dashed curves in the upper panels of 
Figs. | and | There is good agreement with the exact 
results obtained within the sequential-reflection model. 
For TE polarization and the chosen refractive index, the 
analytic results in the literature becomes applicable only 
for Refcr c > 1000. Surprisingly (see however Ref. |l9||), 
in the current situation nice agreement is found by sim- 
ply assuming that the Goos-Hanchen shift is identical 
to the TM case, as is indicated by the dashed curves in 
the lower panels of Figs. || and ^. (For Refcr c > 1000, 
however, it is appropriate to work with the correct TE 
formulas.) 

In conclusion, we investigated the reflection coeffi- 
cients at a curved refractive index boundary by relat- 
ing them to resonances in a circular dielectric disk and 
derived analytic expressions valid in the limit of large 
wavenumbers. The deviations from Fresnel's laws can be 



explained within geometric optics by incorporating the 
Goos-Hanchen effect. In this work we concentrated on 
the wave field confined by reflection inside the microcav- 
ity. The Goos-Hanchen shift also affects the wave field 
outside the cavity [ pl]p^j20[ |, because the modified angle 
of incidence results in a change of the angle of refraction. 
It is desirable to investigate the implications on the re- 
markable emission directionality of non-circular devices 
P],[|, which sometimes departs substantially from what 
is expected from geometric optics pl[ . 
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